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ABSTRACT 

We present the first pilot study of jets from young stars conducted with X- 
shooter, on ESO/VLT. As it offers simultaneous, high quality spectra in the 
range 300-2500 nm X-shooter is uniquely important for spectral diagnostics in 
jet studies. We chose to probe the accretion/ejection mechanisms at low stellar 
masses examining two targets with well resolved continuous jets lying on the 
plane of the sky, ESO-HA 574 in Chamaleon I, and Par-Lup3-4 in Lupus III. 
The mass of the latter is close to the sub-stellar boundary (M^=0.13 M ). 

A large number of emission lines probing regions of different excitation are 
identified, position- velocity diagrams are presented and mass outflow/accretion 
rates are estimated. Comparison between the two objects is striking. ESO- 
HA 574 is a weakly accreting star for which we estimate a mass accretion rate 
of log(M acc ) = —10.8 ± 0.5 (in M yr -1 ), yet it drives a powerful jet with 
M out ~ 1.5 — 2.7 x 10~ 9 M yr -1 . These values can be reconciled with a magneto- 
centrifugal jet acceleration mechanism only assuming that the presence of the 
edge-on disk severely depresses the luminosity of the accretion tracers. In com- 
parison Par-Lup3-4 with stronger mass accretion (log(M acc ) = —9.1 ± 0.4 M Q 
yr _1 ), drives a low excitation jet with about M out ~ 3.2 x 10~ 10 M yr _1 in both 
lobes. Despite the low stellar mass, M ont /M acc for Par-Lup3-4 is at the upper 
limit of the range usually measured for young objects, but still compatible with 
a steady magneto-centrifugal wind scenario if all uncertainties are considered. 

Subject headings: stars: formation — ISM: jets and outflows 
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Introduction 



Collimated outflows are common to accreting objects with a wide variety of spatial 
scales and mass accretion rates, but the phenomenon is best studied in the case of 



jets from the nearb y low-mass stars in their pre- main sequence phase (IRay et al. 



Whelan et al. 



emission lines rtRay et al. 



20071 : 



200 91). These flows present characteri stic spectra rich in shock-excited 



2007 



Hartigan et al. 



19941 ). that allow a detailed diagnostic 



of the kinematic and of the physical properties of the gas. To better understand the 
enigmatic process of jet generation, it is imperative to determine the exact relationship 
between outflows and accretion. Their interplay is at the heart of the widely accepted 
theory of magneto-centrifugal launching , that may be com mon to all jet generating objects 



( iFerreira et al. 



2006 



Pudritz et al. 



2007 



Shang et al. 



20071 ). The pertinent models prescribe 



stringent constraints for the ratio of accreted to ejected mass outflow rates, and it is 
important to understand how this relati onship changes with the mass of the driving sources, 



and in particular for low stellar masses (1 Whelan et al. 



20101). 



The spectra of stellar jets can today be efficiently analysed with dedicated diagnostics 



of key lines, and estimates o f mass outflow rates {M, 



mass accretion rates (M, 



2006 



Hartigan fc Morse 



) (|Bacciotti fc Eisloffel 



out 



1999 



can be compared with corresponding 



Masciadri fc Raga 



2004 



Podio et al. 



20071 ). Reliable estimates of M out / M acc , however, can only be given 



if one has a detailed knowledge of the excitation properties of the gas, and in each location. 
It is thus critical to have a complete spectral mapping along the jet extension, in the largest 
possible wavelength range. 

With these ideas in mind we chose to characterise the jets of ESO-HA 574 and 
Par-Lup3-4, two intriguing objects included in a pilot study of young stars of low and 
sub-stellar mass with X-shooter. Mounted on the Unit Telescope 2 of the Very Large 
Telescope (VLT), X-shooter is an instrument providing high quality, medium resolution, 
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UBV (300-559.5 nm), VIS (559.5-1024 nm) and NIR (1024-2480 nm) spectra. This 
wavelength range includes a large number of emission lines typical of both outflow and 
accretion phenomena, making X-shooter currently unequaled in spectral diagnostic studies. 

The outflows from the selected sources have a continuous jet-like morphology and lie 
close to the plane of the sky, which makes them well resolved spatially with no superposition 
of the blue and red lobes, and with minimum obscuration by the disk. 

Our first target is ESO-HA 574 (llH6 m 03 s 7, -76°24'53"), a low-luminosity, low-mass 



source discovered by 



regi on (d=160 ±1 7 pc 



L w ( Luhman 



Co 


meron et al. 


f 


200 


i) in 


pc, 


Wichmann et al. 


1998) 



(2004 j) in a study of the Chamaleon I star forming 



19981 ). Its bolometric luminosity is only 0.0034 



20071 ). that, compared with other typical T Tauri stars of the same spectral 
type (K8) in Chamaleon I, makes this star underluminous by a factor of about 150. This 
is probably caused by the edge-on viewing geometry of its disk. The lines normally used 



for accretion diagnostic are also very weak (jComeron et al. 



20041 ). which may suggest low 



levels of accretion, yet the star powers a well developed bipolar jet (HH 872, of total length 
of 3150 AU) at positi o n ang le (PA) ~ 45°, revealed in follow-up FORS1 [SII] imaging by 



Comeron fc Reipurthl (120061 ). The jet is made up of 5 distinct knots, with knots A, B, C, D 



forming the blue-shifted jet, and knot E in the red-shifted lobe. 



Nakajima et al 



± 40 pc, 



The second targe t, Par-Lup3-4 (16 fc 08 ro 51'44, -39°05 30.5 ), was first reported by 



(20 00J) as part of a near-infrared survey of the Lupus III dark cloud (d=200 



Comeron et al. 



20031 ). Its mass is only ~ 0.13 M Q , but the accretion tracers in the 
spectrum are quite st rong, and indicate a mass accretion rate of M acc = 1.4 10~ 9 M yr _1 



(IComeron et al. 



2003) . The forbidden lines are associated with a small collim ated jet (HH 



600) at PA ~ 130° ([Fernandez fc Comeron 



2005 



Comeron fc Fernandez 



201ll ). The source 



is under-luminous by a factor of about 25 wit h respect to other M5 young stars in Lupus 



III, its bolometric luminosity being 0.003 L fjMerin et al. 



20081 ). An edge-on disk oriented 
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at about 81° is probably responsible for the subluminosity (IHuelamo et al.ll2010f ) 



We extend these studies, that were limited to the visual band, by taking X-shooter 
spectra of the targets with the slit parallel to the jet axis. We thus obtained in a single 
observation, and for the first time in the field of stellar jet studies, a simultaneous, complete 
spectral mapping from 300 nm to 2.2 /im, that provided immediate and far-reaching 
information about the physical properties of the outflows and of their sources. In this letter 
we present the spectra and discuss the accretion/outflow properties of the two targets. 
A complete mapping of the jet physical conditions derived from the same datasets is the 
subject of a second paper in preparation. 



2. Observations and Data Reduction 



The observations were conducted on April 7th 2010 (Pro gram ID: 085.C-02 38(A)) 



as part of the INAF/GTO program on star forming regions (lAlcala et al 



201lh . The 



seeing varied between 0'.'9 and 1'.'2. In both cases the ll"-long slit was aligned with the jet 
axis, along position angle (P.A.) 44.8 and 130.0 degrees for ESO-HA 574 and Par-Lup3-4, 
respectively (see above). The total exposure time per star was 3600 s (4x900) in all arms 
and nodding observing mode "ABBA" was used (shifting along the jet PA) in order to 
minimize the contribution of the background in the NIR. Slit-widths of l'.'O, 0'.'9 and 0'/9 
were used, yielding resolving powers of 5100, 8800 and 5600 (corresponding to about 58, 34 
and 53 km s _1 velocity resolution) for the UBV, VIS and NIR arms, respectively. The data 
reduc tion was performed ind ependently for each arm using the X-shooter pipeline version 



1.1.0 flModigliani et al 



20101 1 . One-dimensional spectra were then extracted over the star 



continuum and over the extension of the knots. Note that the latters appear only in either 
A or B nodding position, so the spectra over each knot correspond to an integration time 
of 1800 s. The reported radial velocities are systemic, with the LSR velocity of the sources 
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determined using the photospheric absorption line of KI A7698 (VIS arm). For ESO-HA 574 
we measured vlsr = +3 km s -1 , while for Par-Lup 3-4 vlsr — 5.5 km s _1 , in agreement 



with the LSR velocities of the cor responding c 



4 and 7 ±2 km s 1 , respectively ( iDame et al. 



ouds, measured from CO observations to be 



19871 ). On-source spectra were extinction 



corrected assuming for ESO-HA 574 A„=1 .7 mag (derived from Aj =0.45 (iLuhman 



20071) 



using the ex tinction curve at R„=3.1 fromlSavage fc Mathid (Il979l )). and Ai,=3.5 mag for 



Par-Lup3-4 (IHuelamo et al. 



2010 



Comeron et al 



20041 ). Nevertheless, for both targets A t 



values around zero are derived off-source using the [Fell] lines at 1.27, 1.32, 1.64 /im, that 



trace the jet emission (jNisini et al 



20051 ). This suggests that the jet propagates out of the 



envelope obscuring the star, and therefore the spectra extracted for the jet knots were not 
dereddened. 



3. Results 



Figure [T] shows the X-shooter spectrum of ESO-HA 574 extracted at the source 
position, obtained by summing over 5 to 7 pixel rows (0'.'9 to l'.'l, depending on the PSF of 
the considered arm). The system is observed in numerous lines, including both accretion 
and ejection tracers. Balmer lines from H9 to Ha are identified, as well as the calcium 
triplet and He I at 1.083 fim and Pa/3, while Br7 is not detected. Among the many forbidden 



lines, we fine 
As noted by 



strong lines of sulfur 



Comeron et al. 



oxygen, nitrogen and iron in multiple ionisation states. 



(120041 ) . the forbidden [Call] A7321, 7324 lines are much more 
prominent than the calcium triplet, supporting the idea that both the stellar continuum and 
the accretion tracers are occulted by the edge-on dis k, while the jet t racers are unaffected 



(as it occurs in the edge-on system T Cha studied by 



Schisano et al. 



(12009^ . 



In Figure [2] we present Position- Velocity (PV) diagrams of a selection of key lines. 
Knots A, B (jet) and knot E (counter-jet) are indicated in the [SII] A6731 panel. Knot 
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Al, identified in IComeron fc Reipurthl ( 120061 ) as the source, has actually moved to about 
0'.'2 with respect to the star continuum. The emission in Ha and H/3, as well as in 
[OII]AA3726,3729 and [SII] A6731 lines is extended, with higher dispersion close to the star 
and in knot E, that has a bow-shock shape. On the contrary, the emission in the H2 lines 



is con centrated around the star, as often in Classical T Tauri stars (CTTSs) (jBeck et al. 



20081 ). being likely originated in shocked walls of a cavity opened by the flow. In general 
the lines present low systemic radial velocities, with knots A and B blueshifted at ~ -3 
and ~ -6 km s _1 in [SII] A6731 , while knot E is redshifted at +2 5 km s~ l . Tangentia 



veloci ties were determined by comparison with the [SII] images of 



Comeron fc Reipurth 



02OO6I), taken on 21 July 2005. We find V tan = 130, 220 and 325 km s^ 1 for knots A, B 
and E, respectively. Large knot velocity differe nces and velocity asymmetries in op posed 



1994 



Podio et al. 



201 lh . The 



jet lobes are commonly observed in CTTS jets ( iHirth et al. 
determined velocities suggest very low jet inclination angles of only 1.3°, 1.6° and 4° at 
knots A, B and E, respectively, confirming that the d isk is nearly edge-on, and it is causing 



the observed under-luminosity (IComeron et al. 



2004|). 



Turning to Par-Lup3-4, the full spectrum on-source (integrated over the 1" of the 
continuum width) is shown in Figure |3j and three selected PV diagrams in Fig. [2] (bottom). 
Far fewer jet tracers are detected than in ESO-HA 574, while accretion signatures like 
the calcium triplet at 8500 A, are stronger. Pa/3 and Br7 are both present here. Only 
a few lines are extended, notably the [SII] lines, that show a well-developed bipolar 
structure, extending to ~2'.'5 (500 AU) on both sides. Radial velocities up to -20 and 
+20 kms -1 are measured, consistently with the LSR velocities of -18.3 and +23.6 km s _1 



determined by 



Fernandez fc Comeron! (120051 ) in the same line. The strong [OI]A6300 line 



is also extended, although the red-shifted component is affected by subtraction of the [OI] 
atmospheric emission (see inset in Figure [2]). The weak [NII]A6583 shows also traces of 
bipolarity, over similar radial velocities. Tangential velocities have been reported recently 
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by 



Comeron fc Fernandez! (120111 ) to be 170 ± 30 km s 1 in both jet lobes, that combined 



with our radial velocities give a jet inclination angle of about 7°. 



4. Discussion 



A large number of accretion tracers are identified i n the 



Alcala et al. 



spectra at the star position. 



(120111 ). using the relationship 



We estimated mass accretion rates as in 
M acc ~ 1.25 ■ LaccRt/GM* where R t and are th e stellar radius and mass, respectively, 



and L acc the accretion luminosity ( jHartmann 



19981 ) . Th e latter is obtained from the 



luminosities of t he di ff erent lines, via t h e rela tionships in 



Mohanty et al 



(}2005|) 



Muzerolle et al. 



Herczeg fc Hillenbrand 



(120081 ) 



(119981 ). The line luminosities are derived from the 



extinction-corrected absolute fluxes obtained by simple integration in the calibrated spectra. 
Adopted stellar masses and radii, estimated from stellar luminosity and spectral type (see 
Sect. 1) are 0.5 M and 0.2 R for ESO-HA 574 and 0.13 M , and 0.45 R for Par-Lup 3-4. 
The results are reported graphically in the top panels of Fig. HI where the error on each 
point derives from the uncertainty in the calibration relationships and from the noise level 
of the continuum. We considered accretion tracers in all the available wavelength domain, 
including Pa/3 and Br7 (illustrated in middle panels of Fig. H]), while the Balmer jump 
could not be investigated because of the low signal-to-noise in the bluest part of the spectra. 
Excluding [OI], which is dominated by the wind, the determined values for logM acc are 
fairly consistent within the errors. The absence of a wavelength dependence is in line with 



the attenuation caused by an edge-on disk, as the dust extinction 



the disk is not expected to vary in wavelength ( iSchisano et al. 



xom large grains in 



2009b- Excluding [OI], the 



values give an average of — 10.8±0.5 and — 9.1±0.4 (in M yr x ), for ESO-HA 574 and 
Par-Lup 3-4, respectively, where the indetermination is the standar d deviation. Th e M nr , 



value for Par-Lup 3-4 is in good agreement with the one derived by 



Comeron et al 



( 120031 ) 
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The extremely low value for ESO-HA 574, however, can hardly be reconciled with the 
presence of a powerful jet, and the estimate is probably severely affected by the obscuration 
caused by the edge-on disk. Indeed, assuming that both the continuum and the accretion 
tracers are depressed by the same amount, the correction for the disk occultation on L+, 
and, consequently, on R* and L acc , would brin g to a mass accret ion rate of ~1.7 1O~ 8 M 



yr 1 , in line with other stars of similar mass (INatta et al. 



20041 ). On the contrary, for 



Par-Lup3-4 M acc is just abov e the typical va 



ue fo r similar stars, suggesting that the 81° 



20101 ) leaves the accretion funnels visible at 



inclination of the disk plane (IHuelamo et al. 
least in part. 

A comparison of the Ha profile in the two targets is given in Fig. HJ bottom. For 
ESO-HA 574 the majority of the Ha emission is extended and relatively narrow (-100, 
+ 100 km s _1 ), superimposed to a much weaker broad component (-300, +300 km s _1 ), 
localized at the star continuum. On the contrary in Par-Lup3-4 the Ha line is concentrated 
around the star, it is broad (10% width ~ 350 km s -1 ), and its shape is normally not 



associated to a jet component (type II B profile in the classification of iReipurth et al. 



19961 ). In order to disentangle the emission pertinent to accretion or to ejection, the Ha 
line profiles were analysed with spectro-astrometry a technique that allows one to find 
positional displacements of the line peak with respect to the star continuum over angular 
distances of a few tens of milliarcseconds (for details see 



Whelan fa Garcia 



2008|). The 



absence of any significant offset 
jet in Ha in the 2010 images of 



'or Par-Lup 3-4 is cons i stent with the disappearance of the 



Comeron fa Fernandez! (120111 ). and suggests an accretion 



origin for this line. The same applies for the broad Ha component in ESO-HA 574, while 
the jet portion between ± 100 km s _1 has obviously a lar ge offset, preventin g the use of the 



accretion diagnostic based on the width at 10% intensity (jNatta et al 



20041 ) 



Regarding the jets, X-shooter fully reveals its power showing at once the actual 



- 10 - 



excitation state of the gas along the flow. Fig. [2] shows that in ESO-HA 574 strong 
[OIII]A5007 and [SIII] A9531 lines are observed. The knots closer to the source (Al, A) 
present the highest excitation and density as also probed by [SII]A4069, [NII]A6583 , and by 



the [SII] and [NI] lines at 1.03/1.04 /zm. The emission in the helium line at 
very strong, and spatia lly extended, as it has been found in other CTTSs ( 



2002 



Podio et al. 



.083 jLtm is als o 



Takami et al 



20071 ). and 



20081 ). This line is a probe of inner wind regions (IKwan et al. 
indicates high electron temperatures (T e ~ 10 4 K) and densities (n e > 10 6 cm -3 ). In the 
same way in Par-Lup3-4 the absence of lines from high ionisation states, and the strength 
of [OI] and [NI] lines with respect to [Nil] lines reveal that the gas in this jet has a low 
excitation state all along the flow. 

Importantly, this qualitativ e inspection indicates that the usual diagnost ic procedures 



to estimate mass outflow rates ( Bacciotti fc Eisloffel 



1999 



Podio et al. 



20061 ) cannot be 



applied for ESO-HA 574, as the abundances of the species considered (OI, Fell, Nil) 
should be calculated considering high ionisation levels. On the other hand in the Par-Lup 
3-4 jet, instrinsically of low excitation, [Nil] lines are extremely faint and [OI] lines are 
contaminated by the sky lines. Therefore here we give only a preliminar estimate of the 



mass outflow rates (Mout) from the [SII]A6731 
lines are not detected. Following the method in 
calculated as: 

M out (M yr -1 ) = 4.51 x 10 9 (1 + 



ine luminosity, 



Hartigan et al. 



i n the knots where [SIII] 



(119951 ). in each knot M out is 



n c . ^4m(km s *) 

n) W6 731 (L ) Wcm) 



where n c is the critical density for the [SII] ratio (2 10 4 cm -3 ), while Vj an and \t an are the 
tangential velocity and the knot length, respectively. The line ratios and luminosities were 
estimated from ID spectra extracted over the spatial regions marked in Figure 2 ([SII] A6731 
panel), with no correction for extinction (see Sect. 2). The errors are propagated from the 
noise level around the [SII] lines, and from the distance to the star. 
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The results are collected in Table [TJ For ESO-HA 574 n e turns out to be 1750, 
280 and 110 cm -3 for knots A, B and E respectively, well below the critical density. In 
knot Al n e = 3000± 130 cm -3 , but M out could not be estimated because [SIII] lines are 
detected here, while the method assumes sulfur to be all singly ionised. In the other knots 

5-3 10~ 9 Mi? , yr" 1 , within the usual range found for CTTS jets 



M mit turns out 


to be 


(Hartiean et al. 


1995 



Melnikov et al. 



20091 ). Notably, values are similar along the flow, 
as expected if no big sideway losses are present (due, e.g. to large bow-shocks), and, as 
in other asymmetric jets, they a re about the same in the two lobes, despite the gas has 

for a discussion). The obtained M ou t 



different physical conditions (see 



Podio et al 



2011 



however, appears two orders of magnitude higher than the derived M acc , if the latter is 
taken at face value. As discussed above, however, the occultation of the edge-on disk should 
be taken into account. A correction of M acc for the observed under-luminosity would bring 



the mass ejection/accretion ratio back to ~ 0.3 (over the tw o lobes 



for magneto-centrifugal jet launch (0.01 < M out /M acc < 0.5, 



Cabrit 



, in t he expected range 



2009|). 



For Par-Lup 3-4, M out was estimated separately for each lobe, coadding the spectra 
over the two regions shown in Figure [2], bottom panels. The [SII] line ratio provides here 
n e = 1150 and 3200 cm -3 , respectively. The [SII] A6731 luminosit y was evaluated over th e 
same regions, and we take v tan = 170 ± 30 kms -1 in both lobes (jComeron fc Fernandez 



201 lh . Using Eq.dwe obtain M out = 3.6 10^ 10 M yr" 1 for the blue jet and 3.2 10~ 10 
M Q yr _1 for the red lobe, remarkably similar on the two sides of the system. Therefore, 
over the two lobes M out /M acc ~ 0.85, which is slightly above the maximum indicated for 
magneto-centrifugal acceleration. The estimated uncertainties on M acc and M out , however, 
combined with the partial obscuration of the accretion funnels by the close-to-edge-on disk 
can lower this value to 0.3 - 0.5, compatible with that scenario. 



Summarising, X-shooter spectra prove to be fundamental in jet studies as they provide 
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an immediate overview of the outflow kinematics and excitation, as well as a rich accretion 
diagnostics. The first two low-mass (0.5 and 0.13 M ) targets examined proved to have 
accretion/ejection properties compatible with a magneto-centrifugal scenario. An extension 
of the current diagnostic techniques to fully exploit the potential of X-shooter spectra is 
underway and will be presented in a forthcoming paper. 

We dedicate this letter to the memory of Roberto Pallavicini. We thank V. D'Odorico, 
P. Goldoni, A. Modigliani and the ESO staff for support with the observations and the 
X-shooter pipeline, and the referee for very useful comments. E.T. Whelan is supported by 
an IRCSET-Marie Curie Fellowship in Science within the European Community FP7. 
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Fig. 1. — Full extinction corrected (A„ = 1.7) X-shooter spectrum of ESO-HA 574, integrated 
over ~ 1" around the source. 
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Fig. 2. — Position- Velocity (PV) diagrams of selected lines around ESO-HA 574 and Par- 
Lup 3-4. Contours begin at 3 times the r.m.s. noise and increase by a factor of y2. The 
peak value is shown in the upper right corners. ID spectra of the knots in the ESO-HA 574 
jet and the two lobes of the Par-Lup3-4 jets were extracted over the regions marked on the 
[SII]A6731 panels. Distances along the blue-shifted lobe have positive y values. 
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Fig. 3.- 



Wavelength (nm) 

Same as Fig 1, for Par-Lup3-4. The spectrum is extinction-corrected for A„=3.5. 



- 16 - 




Fig. A.— Top: Mass accretion rates from different accretion tracers. Solid line: average 
(excluding [01]); dashed line: ±lcx standard deviation. Middle: Pa/3 and Br7 in ESO-HA 
574 (dashed line) and Par-Lup3-4 (solid line), normalised to their continua. Bottom: spectro- 
astrometric signal in Ha , before (solid curve) and after (diamonds) continuum subtraction. 
Positive offsets are along the blue-shifted jet (dashed lines: spectro-astrometric accuracy). 
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(M yr" 




ESO HA 574 


















Knot A 


1.7 ± 0.1 


7.7 ± 0.5 


130 ± 


30 


3.8 


± 0.4 


1.5 ± 0.4 x 


io- 9 


Knot B 


0.28 ± 0.03 


2.4 ± 0.4 


220 ± 


50 


6.5 


± 0.7 


2.7 ± 1 x 


io- 9 


Knot E 


0.11 ± 0.1 


0.42 ± 0.07 


325 ± 


50 


5.3 


± 0.6 


2.5 i 2 x 


io- 9 


Par-Lup3-4 


















Blue lobe 


1150 ± 20 


1.85 ± 0.3 


170 ± 


30 


7.2 


± 1.5 


3.6 ± 0.6 x 


io- 10 


Red lobe 


3200 ± 50 


5.2 ± 0.5 


170 ± 


30 


9.0 


± 1.8 


3.2 ± 0.3 x 


io- 10 



Table 1: Mass outflow rates in the jets from ESO-HA 574 and Par-Lup3-4, estimated from 
the luminosity of the [SII] A6731 line (see text). 
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